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Abstract 

Contour  integration  is  employed  to  evaluate  analytically 
the  correction  terms  to  the  Rutherford  cross  section  for  high 
energies.   The  sum  of  partial  waves  is  replaced  by  an  equivalent 
integral  which  is  expanded  in  powers  of  the  fine  structure  constant. 
Evaluation  of  the  resultant  integrals  is  performed  wherever  possible 
by  use  of  integral  representations  for  the  transcendental  functions 
appearing  in  the  integrands;  an  unsuccessful  attempt  to  determine  the 
values  of  the  integrals  by  the  saddle  point  approximation  Is  briefly 
sketched.   Those  terms  which  appear  in  the  cross  section  to  the 
fifth  order  in  the  fine  structure  constant  are  explicitly  given  by 
closed  analytic  expressions. 

The  general  form  of  the  correction  term  corresponding  to  an 
arbitrary  power  of  the  fine  structure  constant  is  found  in  terms 
of  two  dimensional  integrals  involving  elementary  transcendental 
functions . 

A  related  problem,  the  non-relativistic  scattering  for  an 
attractive  l/r  potential,  is  also  discussed. 
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1.        Introduction 

h]  hi* 

In  recent  years  several  investigations,  l  J ' L  J  both  theoretical 
and  experimental,  have  been  undertaken  to  determine  the  structure  of 
nuclei  (including  the  proton)  from  the  resxilts  of  the  scattering  of 
fast  (i.e.,  relativistic)  electrons  by  nuclei.  The  experimental  data 
and  the  theoretical  results  measure  essentially  the  deviation  of  the 
actual  scattering  from  that  due  to  a  point  nucleus  and  thus  reveal 
certain  details  of  the  charge  distribution.   It  is  this  scattering 
by  a  point  charge  that  will  concern  us  in  this  paper. 

Since  the  structure  to  be  investigated  has  dimensions  of  the 

-15 
order  of  10   cm  it  is  necessary  to  use  short  wavelength,  i.e., 

high  speed,  electrons  with  the  velocity  (v)  of  the  electrons  close 
to  the  speed  of  light  (c).  Thus,  in  order  to  detennine  the  scatter- 
ing, one  must  use  the  appropriate  solution  of  Dirac's  relativistic 
wave  equation. 

The  appropriate  solutions  of  the  wave  equation  in  terms  of  a 
partial  wave  expension  for  scattering  of  fermions  by  a  static 
Coulcxnb  potential  were  first  derived  by  Mott'-"-'  and  have  been  well 
known  for  masy  years .  The  corresponding  phase  shifts  are  of  a  com- 
plicated enough  form,  however,  that  the  summation  of  the  partial 
waves  in  the  relativistic  case,  in  contrast  to  those  obtained  by 
solution  of  the  non-relativistic  Schrodinger  eqxiation,  cannot  be 
readily  accomplished.  Two  approaches  have  to  be  used  to  determine 
the  differential  cross  section  for  different  energies  and  nuclear 
charges  -  the  first,  essentially  a  numerical  calculation,  and  the 
second,  an  expansion  of  the  wave  function  (or  phase  shift)  in  powers 

_ 

Both  of  these  review  articles  contain  references  to  a  rather  exten- 
sive literature. 
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of  the  two  parameters  Z  (the  nuclear  charge)  times  the  fine  structure 
constant  (e  /he)'-  ^   and  the  ratio  P  a  v/c. 

Bartlett  and  Watson '-  -i  were  the  first  to  carry  out  a  calculation 
of  the  first  sort,  and  their  work  has  been  extended  recently  by 
YadavLJ^  Doggett  and  Spencer  L' -I  and  Sherman'--'.  Yadav  gives  results 
for  uranium  at  energies  of  0.1^5,  001^+,  0.666,  ?.550  and  20.000  Mev; 
Sherman  gives  extensive  tables  for  Z  =  13,  kd   and  80,  for  (3  =  0.2, 
OA,   0.5,  0.6,  0.8  and  0.9  and  for  angles  between  15  and  I65  in 
fifteen  degree  intervals .  Doggett  and  Spencer,  whose  work  was 
approximately  simultaneous  with  that  of  the  last  two  sources,  pub- 
lished resiats  for  Z  =  6,  15,  29,  50,  82  and  92,  at  energies  ranging 
from  moderate  (O.O5  Mev)  to  relativistic  (10  Mev)  over  the  entire 
angular  range.  They  estimated  that  their  error  was  less  tlian  0.5  /o 
except  at  the  largest  angles . 

hfott  L  -J '  L-^-l  carried  out  an  expansion  of  the  phase  shift  to 
determine  the  lowest-order  deviation  of  the  differential  cross 
section  from  the  Rutherford  cross  section  <5p(6)  where: 

?  2 

(A)  6-^(9)     =  ^|-  (l-p^)cscV2 

km   V 

where  m  is  the  mass  of  the  scattered  particle,  v  its  velocity,  and 
e  the  angle  of  scattering.  He  found  the  ratio  R  between  the  differ- 
ential cross  section  and  <5r,(©)  to  be 

(B)  R  =  J\t\     =  1  -  P^sin^e/2  +  nog  sin  e/2  cos^e/2  +  ... 

Following  this,  Sexl  L  -■  using  the  second  Bora  approximation,  foiuad 
instead 

R  =  1  -  P^  sin^e/2  +  nop  sin  e/2  +  ...  ; 
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and  this  result  was  verified  by  the  work  of  Urban L  J  who  used  an 
expansion  of  the  phase  shift .  Dalitz  L  -*  and  Schwinger  L  J  however, 
found  that 


(C)  R  =  1  -  p^sin^e/2  -  «a^  sin  e/2(l-sin  e/2) 


and  this  result  has  been  confirmed  by  the  recent  work  of  McKinley 
and  Feshbach '-  -'  and  Feshbach "-  ''-'  and  Curr  L-  -* .  McKinley  and 

Feshbach  carried  out  an  expansion  of  the  phase  shifts  to  the  fourth 

o 
order  in  a  ;  Z  x  e  /he .  In  their  paper  some  of  the  resulting  series 

of  the  form 

oo 

Zi        f(n)P^(cose) 
0  or  1 

where  evaluated  in  closed  form  and  the  remainder  evaluated  numerically. 
The  most  extensive  calculation  of  this  sort  has  been  carried  out  by 

Q 

Curr  who  expanded  the  phase  shift  to  order  a  and  calculated  the  angle- 
dependent  coefficients  numerically.  By  comparison  with  the  work  of 
Yadav  he  estimates  that  his  error  is  of  the  order  of  3  (a/p)  for 
energies  greater  than  0.5  Mev. 

In  this  paper  we  shall  first  present  an  outline  of  the  Mott 
calculation  and  the  present  a  calcvilation  along  the  lines  indicated 
by  McKinley  and  Feshbach,  and  Curr.  This  calculation  will,  however, 
employ  a  technique  borrowed  from  electrcmagnetic  field  calculations, 
which  is  often  denoted  by  the  generic  term,  the  Watson  transformation"^  -J 
By  use  of  this  last  method  it  is  possible  to  evaluate  in  closed  fonn 
all  the  angular  coefficients  of  the  powers  of  a  that  enter  into  the 
differential  cross  section  to  order  a  in  terms  of  elementary  trans- 
cendental functions  and  the  dilojarithm  of  Euler.  We  then  show  how 
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the  coefficient  of  an  arbitrary  power  of  a  can  be  expressed  in  terms 
of  Sims  of  tvo-dimensional  integrals  vhose  integrands  contain  only 
elementary  transcendental  fvinctions .   Inasmuch  as  certain  of  the  sums 
of  Legendre  fvinctions  times  functions  of  angular  momentum  quantum 
number  contained  in  the  papers  alluded  to  above  converge  very  poorly, 
equivalent  integrals  may  be  more  amenable  to  calculations.   Finally,  we 
indicate  the  nature  of  the  functions  required  to  express  the  angular 
dependence  of  the  higher  order  correction  terms  to  the  cross-section. 


2.   Formulation  of  the  problem 

2.1  The  Mott  formula  L^  -^ 

According  to  Dirac,  the  wave  equation  obeyed  by  electrons  in 
a  force  field  with  a  scalar  potential  V  (vector  potential  A  =  O) 
is  given  by 


(1)      'ir'  ^^^^\   +a.p  +  pmc>  yjf 
where  \^  is  the  four  component  spinor 

(2a) 


In  (l)  p  is  the  momentum,  m  the  mass,  and  E  the  energy  of  the  particle; 
a  and  P  are  the  ^4  x  ^  matrices  indicated  symbolically  by 
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(2b) 


o  6. 

1 


6.     o 


and 


1   0 
0  -1 


Here  6.   are  the  usual  2x2  Pauli  spin  matrices  and  1  is  the  2x2 
unit  matrix. 

Before  considering  the  Coulomb  case  proper,  one  first  develops 

the  formalism  for  potentials  which  decrease  at  infinity  at  least  as 

1  2 

rapidly  as  — —r   ,  e  >  0.   Consider  those  solutions  with  E  >  mc 

r 

that  obey  the  following  boundary  conditions: 


(5a) 

and 

(5b) 


♦-    -finite  at  the  origin 


*. 


— r — >  a.    e 

r->oo       X 


liz 


ikr 


i\(e,0) 


where  r,  9,  (j)  are  the  usual  spherical  coordinates  and  k  is  the 
momentum  divided  byti.   The  differential  cross  section  is  given  by 

(M        d(e,0)  =Z:  |uJe,0)|2y/2:  |aj^ 

Since  it  can  be  shown  that 


(5) 


'5 


the  cross  section  simplifies  to 

i2 

(6)        d(e,0) 


and 


If  the  beam  is  considered  to  be  polarized  along  the  direction  of  pro- 
pagation (z),  then  the  solution  desired  is  of  the  form 


6  - 


^"^^^         *5  7=^^^  ^        ""  ^^^'^^   ^   /"■ 


C^^  %     7=^Sr>  6(e,0)  e^/r 


where  r,  ©,  and  0  are  the  usual  spherical  coordinates. 

ri9i 

By  expanding  the  wave  fimction  in  spherical  harmonics  Darwin'-   -• 
found  the  following  pairs  of  solutions  for  the  partial  waves: 

(8a)  (ij     =  (nfl)P Jco8e)G^(r),        {^J     =  -  G^(r)P^^^(cose)e^'^j 

^  n  O  nil 

(8b)  (%).n-l=  "Vc°«e)G-n-l^^^'      ^^^-n-l"  ^-n-A^^^"""^^^^^^ 


where  G  is  a  solution  of  the  coupled  equations 
n 

dG 

(5,,     i(S.|V.,,,p^,jj^.So^  .  0, 


(5b)      i(f-f--)0„.5^.^F„  =  0, 


after  the  F  are  eliminated.  F  represents  the  radial  wave  functions 
n  n 

associated  with  t-i  and  ^^.  G_  _,  satisfies  the  equation  corresponding 
to  (9a, 9b)  with  n  replaced  by  -n-1.  If  the  asymptotic  form  of  the  G's 
is  taken  to  be  (r  — >  oo), 

(10a)     G„  -  i  sln(kr  -  f- *  \) 

(10b)     G.^_i  ~  i  sin(kr  -  f^  +  ti.^.i) 

then  i°e   G  and  i'^e  "°"  G   ,  are  the  radial  functions  that  will 
n  -n-l 

have  the  desired  asymptotic  form.  Finally,  then,  the  wave  functions 


-7- 


are  given  by 


00   i      i^„     in.n-l 


(lla)    t3=E  i(i>^-l)e  \+ne  ■''■^G_^_^y  i\(cos  e) 


'-n-1 


(m)    t^  =  E  f  e  \  .  e  "°-'g.^.i  iVl^(co8  9) 


i0 


(12a) 


(12b) 


2iTi 
2ikf(e,0)=E  i(i^l)(e   '^-D+nCe 


-'^-i-D   p  (cose) 


^ 


2ikg(9,0)  =  E  (-e  "^  +  e  ^'^"°"^)  P^^^(cos  e)e^^ 
n=l 


For  the  case  of  the  Coulcanb  interaction,  as  in  the  non-relativistic 
formalism,  we  keep  (lla)  and  (lib),  but  the  phase  shift  is  determined 
by  comparison  of  the  asymptotic  form  of  the  G's  with 


i  sin(kr  +  7in  2kr  -  |2  +  ^^) 


where  7  :  Ze  /hv  2a  —  ;  a/p. 

MottL^J    solves  (9)  to  find 


(15a)  e  =  'r-    e 


V^^      Rp+l+i^) 


(15b) 


2iTi 


.      Rp^i+l-i^)      i«(°-P^i) 


Pnfi"^^     Rp^i+i+ir) 


n+1' 


yrp        p  /         2  r— I 

n  -a  ,   r  '  =  7  yl-P  ,   and    |    (x)   is  the  usual  Gamma  function. 


Using  the  abbreviations 


-inp 
ilka)  C     =  -  e         °   r(o„-i7)   /    Rp.  +  1  +  ir) 


(Ub)  F(e)  =  I  H  (-1)°  ^nC^  +  (mDC^A  Pjcos  e) 

(life)  G(e)  =  I  Yl  {-lfL\  .  (nfD^c^^l  Pjcos  e), 

^tott  finds  thatL^^-l 

(15a)  kf(e)     =  -17'  F  +  G 

(15b)  kg(e)     =    Ji7'(l+cos  e)F  +  (l-cose)G    I  /sin  e 

[25] 
where  he  has  used  the  relations  ■-    -• 


(l6a)  P_-,(co8  e)  o  cos  e  P     +  —^  P^-'-^cos  e) 

n+1^  '  n       rH-1       n     ^  ' 


(i6b)  p    -(cos  e)  =  cos  e  p    -  ^i^  p^-'-^(cos  e), 

n-1  n         n         n     ^  ' 


Ejcpanding  F  and  G  in  powers  of  a,   Mott  gets 


(I7a)  F     =     F     +  oF,   + 


(17b)  G     =     G^  +  ojG^  + 


and,  corresponding  to  the  lowest  approximation. 


2  2  2       2 

(18a)     6(e)  =  %^  csc^(|)(l-  ^)(1  -  ^  sinVs) 
4m  V       "     c      c 


and 
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(I8b)     fi-Ce)  =  -^^  csc\|)(l-  4)(1-  ^  sinV  +  m  J  sin  f  ) 

4m  V 


g/v-   2'^    2 
c      c 


corresponding  to  the  next  higher  approximation. 


2.2  The  method  of  contour  integration 


[22],  [25],  [2i.],  [25] 


If  we  have  a  series  of  the  form 


(19) 


S  =  ji  f(n)P^(cos  e) 

n=o 


the  Cauchy  integral  theorem  allows  us  to  write  the  sum  as  either 
.  (    f(v)P  (cose)dv 


(20a) 


2  /    +ijtv 


sin  Jtv 


(20b) 


^  /  f(v)P^(-cose)dv 
2  j     sin  nv 


where  C  is  the  path  indicated  in  Figure  1  in  a  ccxnplex  v  plane  whose 


1/  plane 


Figure  1 

real  axis  coincides  at  the  integers  with  the  values  of  the  orbital 
angular  momentum  quantum  nimibers .  By  the  use  of  this  integral  we 
assimie  that  f(v)  possesses  the  proper  analytic  behavior  so  that 
Cauchy 's  theorem  holds.  One  also  has  the  choice,  of  course,  of  re- 
placing f(v)  in  the  integral  by  any  function  0(v)  such  that  0(n)  =  f(n), 
for  all  integers  n,  provided  the  integral  is  still  defined.  It  may 
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then  be  possible,  by  properly  distorting  or  otherwise  modifying  the 
path  C,   to  evaltiate  the  integral  by  one  of  the  following  methods: 

a)  Saddle  point  approximation 

b)  Sum  of  residues  of  the  fionction  f(v) 

c)  Explicit  eval\iation  by  use  of  integral  representation 
for  f(v)  and  P  . 

The  first  two  methods  are  most  often  used  in  the  study  of  electro- 
magnetic wave  propagation  at  short  wavelength;  the  second  representa- 
tion is  the  resTilt  of  the  so  called  Watson  transformation  proper.  We 
shall  actually  find  that  the  third  method  is  the  most  useful  for 
evaluating  the  relativistic  correction  terms. 

The  most  often  encountered  application  of  these  analytic  techniques 
has  been  in  the  study  of  the  propagation  of  waves  in  spherically  or 
cylindrically  stratified  media.  Specific  calculations  have  been  made 
in  such  topics  as  diffraction  and  reflection  by  conducting  and  di- 
electric spheres'-  -J,  diffraction  by  smooth  objects^  -■,  and  radio 

[2k] 
propagation  •-  -■ .  These  methods  are  most  -useful  in  problems  in  which 

the  wavelength  is  small  compared  to  the  scatterer  of  the  radiation, 
i.e.,  where  k  (the  wave  number  or  2jt /wavelength)  times  a  (the  size 
of  the  scatterer)  is  much  greater  than  unity.  Under  these  circum- 
stances the  simi  of  partial  waves  wl3J.  contain  such  a  larger  number 
(of  the  order  of  ka)  of  significant  terms  that  a  numerical  evaluation 
would  ordinarily  be  called  for;  the  use  of  this  technique  allows  one 
to  determine  the  behavior  of  the  scattered  field  in  an  asymptotic  ex- 
pansion involving  powers  (often  fractional)  of  (r~)'  In  the  mathe- 
matical theory  of  optics  these  ideas  have  proved  invaluable  in  treating 
the  transition  from  wave  to  geometrical  optics  L  -I . 


-  n  - 

In  this  paper  we  shall  apply  the  method  of  contoior  integration 
(or,  alternatively,  the  Watson  transformation)  to  obtain,  insofar  as 
possible,  closed  expressions  for  the  angular  dependence  of  the  rela- 
tivistie  Coulomb  differential  cross  section.  Before  proceeding  with 
a  detailed  discussion  of  the  relativistic  case  we  shall  apply  this 
method  to  the  well-known  non-relativistic  scattering  problem.  Here 
it  is  known  that  the  sum  of  partial  waves  obtained  by  solving 
SchrSdinger's  equation  for  a  particle  in  a  Coxilomb  field  is  exactly 
summable  in  terms  of  a  single  confluent  hypergeometric  function. 
We  shall  see  that  the  method  of  contour  integration,  while  un- 
necessary in  this  instance,  is  capable  of  giving  the  correct  answer 
approximately  without  great  analytic  difficulty.  What  'approximately' 

means  will  be  discussed  shortly. 

[29] 
The  Schrbdinger  equation'-  -^  will  have  the  form 

2         2 

(21)  -^V^t-^t=Eilr,  E>0 

with  the  boimdary  conditions 

\(r  regular  at  the  origin 

, ,         ,       ikr+ir  ^n2kr 

(22)  ^—^^^^^^-^7^Mr-z)^e ^    j 

^      '  r->oo  r  ^    ' 

The  solution  required  iS;,  if  v  again  denotes  the  velocity  of  the  particle, 

(25)      t  »  -^  r(l-l7)e'^/^e"'='  $(i7,l,lk|), 

V  ' 

5  «  r(l-cos  e)  B  2rsin^e/2. 
in  parabolic  coordinates,  and  is 
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(2k)  i  =  -j^   e'^/^  r  n°^-ir)  (2ij,r)^^^  ®n-.l-i7,2r>.2,-21kr)Pjcos  6) 


in  spherical  coordinates.  Here  ^  in  the  confluent  hypergeometric  function 
of  the  first  kind.  The  solutions  here  are  normalized  to  unit  incident 
flux.  By  the  method  of  contour  integration  one  should  be  able  to  get 
fron  (2k)   to  (25). 

Using  the  relations  ^  -* 

(25)  ^(a,c,x)   =  e''^  ^(c-a,c,-x)         JKumraer's  transformation] 

[51] 


and 


1   1    00     1     1 
;r  -  ;r  c  /  ^   -  7:   c-a-  - 


(26)     ^(a,c,x)  =  Hcl.  ,x  ^  2  -  2  M  ^-t^"  2  ^  "  2 ^    L^^   ^jl/2l 


Rc-a) 


in  succession,  one  obtains 


Re  c  >  Re  a  >  0,  Re  x  >  0 


CD 


(27)  t  =  ^z  e^^e'^/-  i:  ^-|—  (2i>fl)Pjcos  e)J2^^(2u)dt 

\/v         n=o  '' 


where  u  =  \/2ikrt .  No  difficxilties  arise  from  the  violation  of  the 
condition  Re  x  >  0  in  (26).  Indeed  the  same  equation  is  used  again 
in  reverse  to  obtain  the  final  answer.  One  can  consider  that 

(28)  2ikr  =  lim   €  +  2ikr. 

€  ->  0 

Actiially  one  need  go  no  further  to  prove  the  equivalence  of  the 

two  relations  (25)  and  (2k).     As  a  special  case  of  a  theorem  due  to 

Bailey  L^  one  has 
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(29)  uJ^(2u  sin  |)     «  i;  {2r^l)J^^^{2u)Tjicos  6); 

using  (26)  with  c  «  1,   a  »  i7,   we  get 

(30)  Mf  -  ^      exp    Q^/2  +  ikz]  Rl-i^)  $(ir,l,2ikr  sinV2). 

\A 

A  theorem  due  to  Erdelyl '-  -•  could  have  been  used  to  circumvent  the 
intermediate  steps. 

It  should  be  interesting  to  see  if  we  can  get  this  result 
approximately  by  means  of  a  saddle  point  evaluation.  Furthermore 
the  discussion  will  seirve  as  a  guide  when  we  consider  the  more 
complex  relativlstic  scattering.  Again,  by  means  of  the  Cauchy 
integral  theorem 

CD 

(31)  F(u,e)  5  JZ   (2n+l)J2^^(2u)Pjco8  e) 

n«:o 

^   [  (2v+l)J2^_^^(2u)P^(cos  e) 

=  2       ; Tuv "^^ 

/      sin  jtv  e 

where  C  is  the  path  of  integration  shown  in  Figure  1. 

Hie  result  of  a  saddle  point  evaluation,  by  virtue  of  the 
assumptions  made  in  the  method,  holds,  strictly,  only  asymptotically; 
that  is,  if  f  (u,0)  is  the  result  of  a  saddle  point  approximation  for 
the  fvmction  f(u,e),  f(u,e)  -  f  (u,e)  _   >  0.   In  other  words,  the 
first  term  in  the  asymptotic  expansion  of  the  integral  in  (31)  (as 
obtained  by  the  method  of  steepest  descent)  should  be  equal  to  the 
first  term  in  the  asymptotic  expansion  (large  u)  of  uJ  (u  sin  6/2). 

We  shall  assxane  first  that  a  saddle  point  lies  in  the  first  quad- 
rants and  that  Im  v   »  0,  where  v   denotes  the  position  of  this  point, 
sp  sp 


Then 
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(52a)  e       sin  rtv  j^^_^^  >     -     ji 


(52b)  P  (cos  e) 


Inn^>oo 


->^/=S=4-i'-l 


1  .  1 


(3.C,    n-i.  .  ^}:iiul  ,^>  (. .  i,-/^ 


Rv+i)     r(v+|  +  i) 


2 


and 


35] 


2iusinX-i(2v+l)X  ^^       2V+1 

(32d)  J,^^^(2u)  ->  ^ ,  '   -ost.^-^. 

Using  these  relations  in  (51), 

(  iR  iR 

(55)  F(u,e)     =  ^  \/cosXsln  X     (e     "^  +  e     ')d  X 

c'    vn  i  sine 


where 

R^  =  2u  sin  X  -  (2v+l)  X±(v  +  |)e+f 

=  2u  sin  X  -  X  cos  X  ±  2  ^°^  ^  I  ^  f     • 

The  position  of  the  saddle  point  in  the  first  quadrant  is  given  by 
the  solution  of 
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1       s 
and  is  X  gp  =  ^/^^  ^gp  +  2  ^  ^  ^^^  2   '     '^®  equation  of  the  fall  line 

or  path  of  steepest  descent,  is  given  by 

(55a)   Impiusin  X+  i(2v+l)  1+   i(v+  |)e  -  ^=   Im[2iusin  X  sp+  Q 


or 


(55b)     Rej^sinX-  (v+  |)(  X-  |^  =\/krt  sin^  e/2  . 

In  the  neighborhood  of  v   we  place  {\^\      +  ds  e  ^)   and  obtain 


Re  +  e  '    =  0. 

L   v/p         J 


The  solution  0  =  it/8  yields  an  exponent  which  decreases  as  one  moves 
away  from  the  saddle  points.  The  results  for  R_  are  obtained  in  a 
similar  fashion. 

The  contribution  to  f(u, 6)  from  the  saddle  point  evaluation  is 
given  by 

cos(2usin  -^  -  r) 
(36)         f   =  u  —      ^   ^—  . 
^^         v/rtusin  e/2 

Since 

(57)  J  (2usine/2)  -5—- >         ^  cos(2usine/2-nA)+   0  (-^       ) 

Vjmsine/2  u'^/'^ 

we  have  f(u,e)  -  f  (u.e)  — r >  0.  We  shall  employ  the  method  of 

sp      u— -Xd 

steepest  descents  in  order  to  estimate  the  relativistic  correction 
terms  in  the  next  section,  but,  unfortunately,  the  results  are  not 
useful . 
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2.3  Relatlvlstic  wave  fiinction  and  method  of  steepest  descent 

Here  we  consider  not  the  phase  shifts  but  the  wave  functions 
themselves.  Mott  finds  for  the  G's  : 


(38a)     G.^_^  «  I  ^_:r^ (2kr) 


^^-ikr^Ky/2  p^    iRPn+l+i^)! 


rr(2o^+l)        v/Pn-l>   f-iVn-i^' 


X  <(Dj^-i7)0(p^+i7,2p^+l)2ikr)-(n-i7')0(pn+l+i7,2p^+l,2ikr) 


and 


T   -ikr  Tcr/2  o-    |r(o-+l+i7)l 

le    e'   /^,  xn   ''^n     '' 


(38b)     G^  =  I  ^^r-^ (2kr) 


r  n(  20^^+1 )        W^^  Vn+l+i7  ' 


X  <(o'^-i7)0(cJj^+i7, 20^+1, 2ikr)+(n+l+i7')0(<Jj^+l+i7, 20-^+1, 2ikr) 


where  ff  «s  p  ,  . 
n    n+1 

It  is  convenient  to  deviate  from  the  procedure  of  previous  authors 

by  splitting  off  the  term  that  in  the  limit  c  ->  oo  yields  the  non- 

relativistic  wave  function.  For  this  purpose  we  use  the  identities'^  -■ 

(39a)   0(p^+i7,2p^+l,2ikr)  «  0(p^+l+i7,2p^+l,2ikr)  -  |^^  0(p^+l+i7,2p  +2,2ikr) 

n 

and 

a  +i7 
(39b)   0(a^+i7, 20-^+1, 2ikr)  =  _°  ^^^  0(o^+l+i7,2<y^+l,2ikr) 


20^ 
^  ^5^  0(Vi'''2^n'21i^^) 


In  addition  to  this  we  use  an  integral  representation  for  the  0's 
mentioned  previously  in  this  paper.   One  finds  that: 
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00 

X  J^p    [2u]dt  +k  exp[%l  -inp^+ikrf  g|  e'V^^J^p  ^ifSuldt 
°  o 

-  k  exp[iki^  f  .  if  -ino-J J     ^^—  J^^  .^(2u)dt. 


In  the  non-relativistic  limit  (y'   ->  7,  P  ->  n,  etc.)  the  first  term 
in  both  (UOa)  and  (40b)  approaches  zero. 

Substituting  these  expressions  into  (lla)  one  finds  (exchanging 
the  order  of  summation  and  integration). 


/^  c   -t  -i7  00  -Irtp  -Incf   I 

(.1)   ,^  .  .e---/^  /  ^-^  E  kp^,,(2u)e    ^-(..l).,,^.,e    ^j 


(-i)^^(co8  e) 

.   icy/2  +ikr  /    +  <^  ,  od     ^       \  r-  — i 

+  1  ^-4 I   e-V^^-^  i:(-l)\(cose)/[^n.l-c.^)+i(7'-7^ 

-irtC  p  ^   "^"Pn     I 

Xe    ^J2^(nfl)+|(^p^-n)+i(7'-72jne    ""j^^    \     . 

We  next  expand  the  Bessel  functions  of  non- integral  order  about  the 
integers  to  obtain 
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(42)   ^3  =  l.e^^-^^^/2  [  ^-f^  E  L  ^&P,-n A(n.l)(cr-n.l)^ 


n=o  nt=o 
o 


on 


^  ^ilaM.:r7/2  /°°.t  -i^-i  0°  ^  1 


2     r 

'o 


n=o  m=o    ' —  — '  ' —  — 


^  iJ^e''^hg-n^l.iiy^-y)\^^-n^T\^ 


3n' 


^"^        -irtn, 


2—  (J^^e"^'^)  W-D'V  (cos  e) 
v^m  '^  2IW-2     '  [  ''   '  n^ 


We  shall  now  see  that  the  relativistic  correction  terms  cannot 
be  obtained  from  the  saddle  point  contribution.  An  inspection  of 
the  terms  indicates  that  the  corrections  arise  principally  from  the 
partial  waves  of  low  angular  momentum  or,  in  mathematical  terms,  from 


\/2~~2 
/n  -a 


those  terms  in  which  /n  -a  differs  significantly  frcxn  n.  The  saddle 
point  previously  obtained  in  dealing  with  the  non-relativistic 
solutions  corresponded,  on  the  other  hand,  to  large  ccanplex  i/alues  of 
V,  and,  strictly  speaking,  the  use  of  the  method  of  steepest  descent 
assumes  that  u  and,  hence?  v  are  large  . 

The   use  of  the  method  as  presented  here  is  of  doubtful  validity  inasmuch 
as  an  integration  over  all  values  of  t  from  zero  to  infinity  must  be  per- 


formed ss  indicated  in  (2?)  and  u  =  /2ikrt.   The  contribution  from  the 
interval  0  to  l/2ikr  cannot  be  ignored.  Asymptotic  expansions  of  0(a,c,  x  ) 
for  a,  c  and  x   all  large  would  be  useful  in  this  connection  as  the  integral 
representation  would  not  have  to  be  introduced.  However,  these  are  not  known 
to  the  writer .  [Those  quoted  in  HTF,  Vol.  1,  p.  28l-2  are  not  applicable.] 
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We  shall  consider  only  one  term  of  the  wave  function,  nwinely 


00  ^^'.-t  -ir 


(i.3)   h(r,e,a,7)  5  ke^^^'^'^/^  JZ     j  ^—^-   [^(p^-nK(nfl)(o;-n-lT] 


(-l)\(cose)  ^  \£'^   J2„.i(2u)]  dt. 


If  an  expansion  of  the  radicals  p  and  cr  is  made  in  powers  of  a  , 
the  first  term  would  be 


(44)   -  |-ke """'/^ 


»   '"e-V^^ 


2E    S_i-(-l)\(cos94 


•inn 


GO  .   . 


-"""^2«l(^") 


dt 


Using  the  asymptotic  expansion  (32b)  for  the  Bessel  function,  one  has 
(lis) 


Sj(»"'""^2v.l'~E=^X-l3 


-ijtv- 

e    J, 


2v+l' 


and  repeating  the  procedure  used  above  for  the  non-relativistic  case, 
the  method  of  steepest  descent  yields  (considering  a  and  y   to  be  in- 
dependent ) 


(kS)       h-^^)(r,e,0.v).^h-(r.e.a,r^    ^  1.  exp  rikr^fllAcosf     ^ 


Ba' 


-t.-i7 


V'misine/2 


X    |(-ie-irt)exp[2iusin|  -  J^J  +  (ie-i«)exp  [-2iusin|njj 


dt. 
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The  prine  indicates  that  the  last  term  in  {hk)   has  been  drorped. 
These  integrals  can  be  performed  in  terras  of  parabolic  cylinder 
f\mctions;  one  can  then  find  the  behavior  for  large  kr.  Using  •-  ■^'  '-•^-^-' 

(^7)   j   e-^^'-l  exp[I(2at)l/f!  dt  =  e^/^  2^^'  ri^v)^,^/.^^   , 
o 

Re  v'>  0  , 


fM 


1         -l/ke 


(i+8a)  D  .U^^^)^—>     i-  e 
'   -2v^  'a->oo    v 


Ij—  <  2  arg  a  <  -  5Jt/4  , 


1/2    -2iTCv' 


(hf^)     n   r^l/^'i ^s  1   .-lA^    (2rt)^^-   e"'^^"^     rt^l       ^^ 

(48b)  D_2^,(a   )^-^>  -  e       '  ^T^-      -^T/T       ,  I^  <  2  ^^6  ^  <  t"  ^ 


(2v')    a^ 


one  has: 


/-,  \             kexp  [ikr+rtr /2] 
(49)      h'^-'\'r,G,0,r) 


^  V''it  cos-  \/2ikrsin  ©/2 


x<^ 


(-IG-iit) 


expLiikrsir.^e/2]   2-^'      "*"   ^^    p  (5/2  -  2iy) 


(ie-in)2^A  +   i^exp  [^ikrsin^  rt/g]    f  (3/2  -  217) 

(ie-irt)exi,r2ii"sin^e^    (2n)^/^  2^/^  ■"  ^^  2^(7-  Ji) 

_I IZ e 

(^e       '    krsin  9/2)    ' 


TliiE  js  clearly  incorrect;   first,    it  is  not  of  the  foiTi    exp[i.kr/r]; 
secondly,    it  predicts  an  infinite  cross  section  in  the  back  direction  (9=rt) 
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In  point  of  fact  this  term  shoixld  yield  the  terra 

nai3  sin  e/2  (l  -   sin  e/?) 
in  equation  (C)    and  it  will  be  shovm  later  that 


Jl).      ar.n^ ^.       exp  [jkr+iyln2k£j  2     ^_^^^  ^^^ 

(50)      h^  '{r,9,0,C)-^^> T sine/2 


For  these  reasons  we  will  adopt  a  procedvire  equivalent  to  that  of 

2 

Mott  and  others  and  expand  the  phase  shift  in  powers  of  a  and  y .  Having 

done  this  we  shall  use  the  method  of  contour  integration  to  sum  the  re- 
sulting series . 

3.   Calculation  of  cross  section. 

5.1   Expansion  of  the  scattered  field  in  powers  of  the  fine 

structure  constant. 

For  convenience  each  component  of  the  wave  is  spljt  into  three 

parts.  The  expressions  for  the  direct  wave  are: 

oo 


(51a)  J^)  =  e^^-  e«^/^  (    ^1^   Z   ( -1)^  :^(cose) 

c 

X       nexp[-ijtpj  J       ^^  -  (n+1)  Jg^  .^  exp[-irtO-J   dt 
'—  a  n  — ' 

(51,)     4=)  =- iZliZi    e^      e-^/S     [       ,-\-iy-l    f;    (.D-rtcose) 


n=c 
o 


X       exp[-inpj   nJ^^     +  (n+l)  Jg^  exp[-irttfjj      dt 
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X 


n(o^-n)  Jg^^       exp[-ijtpj    -   (n+l)  (cJ^-(n+l)  )J2^    expL-inCTj 


dt 


These  equations  are  obtained  directly  from  {^l)   but  in  doing  so  we 
have  dropped  a  factor  of  k  so  as  to  later  obtain  an  incident  wave 
of  absolute-  value  unity.   The  results  for  \tf,  can  be  written  down  by 
inspection  of  (lla),  (lib)  and  (31).  One  sees  that  the  changes  to 
be  made  are  as  follows: 

a)  replace  F  (cos©)  by  P^  '(cosS)  and  multiply  by  e  '^  ; 

b)  drop  the  factors  of  n  and  n+1,  multiplying  the  B'3ssel  functions; 

c)  reverse  the  sign  of  the  factor  containing  J-     or  J      • 

n        n 

The  use  of  three  separate  wave  functions  is  equivalent  to  writing 
(l3a)  and  (ijb)  as 


(52a)  exp[2i^_^_£|  = 


r(o„+i  -  ±y) 
r(p„+i  +  ±y) 


exp[ijt(n-P^)] 


l+i21II 


p  -i7   P-i-y 
n      n 


(52b)  exp[2i^J  = 


RV^-ir) 
exF[ln(n-cf  )] 


1+i 


n+1-0-  - 

7-7 '    n 

o-^-i7  "^  <5. -i7 
n      n 


The  Bessel  functions  in  equations  (51a,b,c)  and  in  the  corresponding 
for  t^  are  now  expanded  in  a  powers 
we  obtain,  for  example  [see  equation  U2] 


ones  for  t,  are  now  expanded  in  a  powers  series  about  the  integers  and 
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,^,,   ,(1)    ikr  n7/2  e  t     <;—  't-  n 

(55)  t^  '  =  e    e  '    I 1_  1.        —i 

'  n=o  m=o 

o 


X 


X   [Pn-^-l]' 


2 
We  consider  that  the  factors  such  as  fp  -n]  are  expanded  in  powers  of  a  . 

The  t  integration  is  carried  out  in  all  terms  except  the  coefficient  of  the 

2 
lovfest  pov/er  of  c,  in  each  of  the  three  parts  of  \I',  and  of  ilf.  .  Taking 

the  asymptotic  form  for  large  kr  of  the  resulting  expressions  yields,  with 

the  exception  of  the  six  terms  [the  first  three  terms  of  -^^^   and  ^I/H  mentioned 

above,  the  form  of  the  incident  and  scattered  wave  far  frcxa  the  nucleus. 

This  procedui'e  is  of  course  entirely  equivalent  to  a  phase  shift  expansion 

and  as  such  is  really  the  same  method  as  that  employed  by  previous  authors'-  ^'^    ^ 

with  this  exception.  The  first  terms  in  -iiX    '     and  \}r^  '  contain  sums  of  Bessel 

functions  that  are  exactly  sumnable.  No  futher  expansion  of  these  terms 

(in  powers  of  y)  is  necessary.   The  first  term  of  an  expansion  of  these 

terms  in  powers  of  r  agrees  with  that  of  McKinley  and  Feschbach  •-  -' .  ( Cf . 

their  equation  (6)ff.  with  our  equation  (56c)  ). 

We  shall  concentrate  first  on  these  six  terms  involved  and  denote  them 

in  generally  \!f ' .  One  gets 


-   2Jf 


.-t.-ir 


(5i.a)  ^l^^^     =  e^^e'^/2     |       £X_      ^  (2n+l)J2^^[2u]   PjcosG)  dt 


=  e^^e'^/^     j       ^-t^-ir  j^(2u6ine/2)  dt 


ikr 


/ci.-u\      >  '(2)  1  /         ,\  JO'/2  e 

(54b)     t^^        =  ■  2  ^'^"^^^  ~kr 


00  ,  V 

o  I  J 


X       P  (cose)  e"^"^''"-'-  dt 


—   i-/-t'i    p      '       — s-In— 


I  (r-7')  e^^/'^-S^^     sin|       J       uJ^(2usine/2)e""t 


t,-ir-i 


dt 


o  L.-'-  - 

CO 

=  -  1^     e'^/2  ^2     I       e-S-i7-l    [jj2usine/2)   -  Jj2u^   dt 


dt 


(5^d)      irj^^l)      =  0 


Ikr 


!(,.,.)  ,«r/2e^  ,10 


j     ucGse/2  J^(2usine/2)e"H"^^'-'-dt 


=     ctne/2  e^^  i^^^ 
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^5^f)      h^^^   =     k    exp&  r/2]   a^     [        e-^-i^-1 


CO   J,  P^l) 
r—     2n  n 

1  ^ 


00       J 


-  i: 

o 

i(ZS     i 


2n+l  p(l) 
n+1        n 


-t, -iv-1 
e      t  dt 


2  f       -t^-iv-l 


2kr 


tane/2exp[Tt  r/2]   a  e"   t' 

o 


J  (2usine/2)-J  (2u) 
o       '  '   o 


dt 


=  -tane/2  rj'^   e^^ 


In  order  to  obtain  the  asymptotic  behavior  of  these  terms  for  large 

[40] 
u  we  use  '-  -^  . 

oo 

(55)  [   e"Vj  .[2(xt)^/£]dt  -— — •>  ^-^v"-^^vrQ^':K/2^->  r  (P/g+i+g)  (x)-oc'-l 

i       1         2 
and  thus  obtain,  denoting  the  asymptotic  form  of  t  .  by  u  .,  In  sin  9/2  by  X 

and  P  il-iy)/P   (l+ir)  by  M(l) 


(56a)      u^"*-^   =  exp[i^kz-i-/ln2kr  sin  e/gj    + 


2ksin  e/2 


M(l) 


.    ,    exp  rikr+ivln2krl 
ivX       ^  ' —  — ' 


,/^x  ,      ■    y    exprikr+iyln2krl 

(56b)   u;^2)  ^  .  ^(^)  2:^  ^xn  _^L_ J 
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,  .       2          ±y\      ,    exp  rT-ylngkr+ikrl 
(56c)  u;^5)^    a_   ^^^)   ^_:1    __L_ -J_ 


(56d)  uj^^^)  =  0 

(56e)  u;^2)  =  ctn  e/2  e^!^  u;(2) 

(56f)   u;^5)  ^  .  ^^^e/2  e^i^u^^^) 

exp  jTkr+i7ln2kr] 

The  coefficients  of  in  (56a,b,c,e,f )  shall  be  denoted  by 

r 

'    '   2  '   '  2  ' 
f  ,  fp,  a  f   ,g^,(x   g  respectively   In  applying  (55)  to  obtain  the  behavior 

for  large  kr  for  ilr^^'^and  \lf|  ^-^  ,  one  should  consider  J^(2u  sin  ©/2)  -  J  (2u) 

^  H-  o  o 


as  Lim   J   (2u  sin  6/2)  -  J   (2u)  in  order  to  satisfy  the  validity 
€->o   ^  ^ 

conditions  specified  by  HTF. 

In  order  to  evaluate  the  remaining  terms  in  the  scattered  waves, 
one  applies  (55)  to  (^2),  term  by  term,  and  in  this  manner  obtains 
series  of  Legendre  functions  with  coefficients  involving 

P  (n+l-iy)/!""  (n+l+lr).   This  fionction  is  expanded  in  an  ascending 
power  series  in  r .   For  v  ~  c,  7  is  approximately  equal  to  a  and  therefore 
terms  of  order  a^7^,  say,  must  be  lumped  with  terms  of  order  a  .   If  we 
denote  the  resulting  series  of  Legendre  functions  occurring  in  u^  as  S. 
and  those  in  u,  as  T.  we  have,  for  example 
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I  2        2 

(57)  2ku^  =  J  2kf^  +  2kf2  +  2ka^f^  +  ^  ^1  +  ¥  Pj  "  ^^'^^gl 


-  -iZllJ.  [S3  .  i^S^  -  -^   [T  2«S,  .  lilS^ 


^LiZlll)  ,[2iS^^  -  2iS^^  -  .Sg  +  2«S^3  .   2iS^^ 


+  g-   [5«Sg+  ln\-  213^7]+...  1 


exp  [Ikr+i'>'ln2kr] 


The  result  for  u,  is  analogous  with  the  f s  replaced  by  the  g's  and 

the  S. 's  by  the  T. 's.   The  T. 's  can  be  written  down  by  inspection  from  the 
111 

S  's  by  a  method  equivalent  to  that  by  which  i|/>  was  obtained  from  >|r  .   The 
series  S.  and  their  sums  as  obtained  by  contour  integration  are  tabulated 
in  Table  1 

The  meaning  of  the  symbols  in  Table  1  is  as  follows: 

(58)  *^  (n+1)  =^^     log  r  (n+1)  > 

dn 

,z 

/   =  -     ln(l-|)d|/^,   the  dilogarithm  of  Euler  , 


1-sin  9^2 
^  ~       1+sin  e/2 


3 -2   The  evaluation  of  the  series. 

As  an  example  of  the  evaluation  of  the  S  s  by  contour  integration 
consider  S,. 
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(60)        S     =  ^       tgCn+l)   Pjcos   e)   +  ^     tgCn+l)   Pjcos   9) 


^  t2(n+l)Pn^cos  0)  -  t^^^) 


(61)   XI  ia^''^^^  Pn^'^°^®^  =  I 


p  (-cos  e) 

: ■i|r_(v+l)dv 

smrtv       2 


O^ 


,     -  ^„(l/2  +  iT)dT 
coshrtT      2'  ' 


In  this  last  step  ve  have  distorted  the  path  C,    as  indicated  in 
Figure  2.  so  that  the  path  of  integration  lies  parallel  to  the  imaginary 
v-axis  after  adding  the  circular  arcs  at  infinity. 


\ 


r 


V  plane 


/ 


/ 


Figure  2 
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TABLE  I 

Part   I  -   Sims  appearing  in  cross-section. 


(P    =  P  (cos  e), 
n         n 


CD  00 

S,    =  y~     P     +  y~  P     =  CSC   e/2   -   1 
1       •^       n       '^—     n  ' 


Sp  =  ^  P     t-,(n+l)   +  ^     P     ijf-,(n+l)   =  \|/(l)  [csc9/2-l]-csc   e/2[jn(l+csc   e/2)-lnt] 


OD  00  r~2  -/  -/ 

-^1  n  I —  — 


-  //6 


S,    =  ^  P     i^(n+l)   +  Z!  ^n  ^^^°+^)   =  ^^^  ®/^    [i(l)+«n(i+csc   G/2)-ini+]^ 
1       "  o  ' — 

S     -?     !2  +  y-     !2 2^      / ^^      -U.-^        /2sine/2\      _^£      /l+sinGM 

^1       n^        o        (n+lf     '         ^    il+^i^®/2)  2      il+sine/2J        ^2^2         j 


-   6^^(l-s±n  e/2)    -   en2(l+sin  e/2)   +  ^in2in(l+sin  e/2)    -  2in   2 


2  ^TT 

-  2in    sin     e/2  Pm{1+  sin  e/2)   +  ^-     +  S. 


OD       P  oo      P 


_      _         1        S  JtV 

2         ^     ,      ,  .2     ~  "  2       .2  sinrtv   ^v 


In  o        (n+1) 


7T7T7::T:P,(-cose) 

ov 


/^(cosVs)    -  f- 


00  P  CD        P 


(1+sin  e/2) 


9       ^       n  ^  (n+1) 

oo        P    l|f,  ( n+1 )  00        P    l|f,  ( n+1 )  J'  /r^        •        r,/r,      \  -P 

^15  -  2_  n  "  2_         n+1  -     ^     2      \^1  +  sin  e/2J      "^     ^\L^\rMjd) 

-  c^p(m)    -  ^„(-m)    -   2i(l)£n(l+sin  9/2)+  inUin  ^i  +   sin  e/2) 


£n     ( 1+sin  e/2)+   2^n 


i-iilLe/Liin 

1+sln  e/2' 


1+sin  9/2  "~1  ,  7rt^ 


U-sine/2)^'  ^ 


>, 

Ti 
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H 
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eNo 


sUl. 


H|CM 
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v^ 


8U|o 
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^ 


CD 

fi 

•H 
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W 
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CM 

CM 

H 

C 

v.^^ 

"S^ 

G 

^ 
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1 

,_^ 

«|vo 
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V 
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"^^io 


^ 
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H 

eNo 
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8Nr 
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Part  III  -  Stuns  appearing  in  [i^^. 


In  general,    if    ^  .   =  T.     ^^"^^  PjcosG)  +  ^^  g(n)  P^(cos©), 
1  o 

thenT^=f:     fMp^(cos©)?     X     ^^n^^^^®) 


T     =  -  tan  9/2  S^ 

T-  =  -  tan  e/2  S^  +  cscG  S^ 
T  =  -  tan  e/2  S  -  esc©  Sg 
T,  =  -  tan  e/2  S^^  -  esc©  Sg  +  2csce  S^^ 

T  =  -  tan  ©/2  S  +  csc©(S^  -  S^^  +  S^) 

Tg  =  -  tan  0/2  S^ 
T  =  -  tan  ©/2  S^ 

To  =  ctn  e/2  Sg 
T^  =  ctn  ©/2  S^ 
^10  =  ctn  ©/2  S^Q  -  CSC  ©  S^^ 

\l  =  -   tan  ©/2  S^^ 

T^  =  Ctn  ©/2  S^  +  CSC  ©  S^^ 

T   =  ctn  ©/2  S   -  CSC©  Sg 
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The  "behavior  uf  the   integrand  on  these  circular  arcs   is  dominated  by 

p  (-cose) 


sin  rtv    I  Im  V  I  ->aD 


->  exp 


-©  Im  V 


for  0  <  9  <  n  and  thus  the 


integration  over  the  circular  arcs  yields  a  vanishing  small  contribution 


to  the  integral.   Use  is  made  of  the  following  integral  representations 

-°°        cosh|Xv+l/2)t]dt 


M'M 


(62a)    P  (-x;  = 


r  (1/2) 


^    nv+i)r(-v) 


\/cosh  t  -  X 


Re(v)  <  0,  Re(v+l)  >  0,  +  x  not  on  the  real  axis  between  1  and  en  ; 


-t  -tz 
(62b)   ^^^(z)  =t^(l)  +  I    ^  "f^ dt  ,  Re(z)  >  0  ; 


,00     -tz 
(62c)  t^(z)  =   I    i^  dt,  Re(z)  >  0  ; 
1-e 


in  evaluation  of  the  sums, 


Since 


:63) 


costT  dt 


\/cosh  t-cos   0  I         Vcosht-cos   9 

O  -QD 


e      dt 


we  can  write  the  sum,  using  the  integral  representations  as 


5^ 


(6k)  J2       ^2^"""^^^    ^JoobQ)   =  ^ 


2        I  I  I  |e  ^    e         e 


-OD         -oo        o 


( l-e~    )  \/cosli  t-cos  9 


d|dtdT, 


Interchanging  the  order  of  integration  and  applying  the  Fourier  integral 
theorem,  we  have 

-S/2 


o  ^  )       Vcosh  I -cos  e  (1-e  ^) 


^1"W losing  X  =  e-g 


o    [(x-cose)^+sin^eJ^/^(x-l) 
G-rt/2 


;e/2 


In 


cos( v-6) 


dv 

sin  v-e/2 


e/2 


using  tan  v  = 


x-cos( 


sinG 


Hence,  by  use  of  the  further  substitutions  u  =  v-e/2,  t  -  tan  u/2 


tan 


(66)      2  i:     ^r^Cn+l)    PjcosG)    =  +   esc   G/2      j  ,n  {|^^}f|^|   f 


nr 
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where  c  =  tan  9^2  +  sec  ©^2,  d  =  tan  G/2  -  sec  ©/2  . 
Finally 

■  —    2 


(6?)   S  =  CSC  e/2 


j£  _  /      A-smeM    +  /,    /    i-sine/2 

~     2   ll+slne/2J     2  1  l+sine/2 


For  some  further  details  see  Appendix  1 
Aq  Interesting  example  Is 


00  P  (cose)    00   P  (cose) 

(68)  s.  =  i:  ^L^. —  Y.    " 


2      ^-—       2 

In       o    (n+l) 


r 


■^  2   J    slnrtv  ~   v^(v+i)2  ^^   '   ^^     )  3     slnTtv  '^^ ' 


I   Is  here  a  circular  path  about  the  origin  of  radius  less  than  l/2, 
say.  taken  in  a  clockwise  direction   When  the  integral  over C Is 
distorted  to  run  parallel  to  the  imaginary  v  axis,  it  will  have  the 
factor 


T  cosh  ItT  dT 


=  0. 


Evaluating  the  residue  at  zero  of  the  Integral  along  [~*  one  has 


(69)   S.  =  -  i  ^    f'p  (-cos  9)-^       I      =  t(oos^Q/2)   -  i 


2    2    \,  V       'smrtv  j       \       _  2'     '  -^   6 

3 -3  The  cross  section. 

Since  the  leading  term  in  the  scattered  wave  amplitude  Is 

proportional  to  a(for  a  ~  ■^),  knowledge  of  the  wave  function  to  terms 

k 
in  a  permits  the  calculation  of  the  differential  cross  section  up  to 
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5 
those  terms  proportional  to  a  .   However,  we  are  actually  interested  in 

the  ratio  of  this  cross  section  to  the  Rutherford,  cross  section, 


2  2  2 

(70)  Op(e)  =    p"'  i, =  ^ol      (l-3^)csc  e/2  [Relativisti£| 

^       4k2sin^e/2     l^mV 


The  ratio  c(e)/o^,  commonly  denoted  by  R,  will  contain  terms  proportional 
to  a    if  the  wave  function  and  cross  section  are  known  to  order  a   and 
a    respectively.   We  shall  indicate  the  order  in  a  to  which  the  wave 
function  has  been  calculated  by  a  Roman  numeral  subscript.   Then 


(71a)   R  =  l-p^sin^e/2 


(71h)   ^ii-^i  =  «°-P  ^i   sin^e/2  =  nop  sin  e/2  (l-sin  e/2) 

2  2  2  2  2 

(71c)   Riji-Rji  =  ^4"  sin^e/2  +  ^  ^  ^   S^  sin^e/2  sec^e/2 

2  2 

+  ^^j^     \^sin  e/2   sec^e/2  +  a^sin^e/2(2S  +P^Sq) 

(71d)        Rjv-Riii  =  ^     q'>^S^   sec^e/2   sin^e/2   -  ^  Q^S^   sin^e/2 
+  rt  ^     S^(2S  +Sg)sec^e/2   sin  e/2 

■Z.  T 

+   rt  ^  (S^-2S2)\sin^e/2   sec^e/2  +  ^    pS^-p^sH  sin^e/2 


+  2rtx^pS^      sin^e/2   -  ^—  S^sin^e/2   -  ^^     sin^e/2   Sg. 
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Here  Q  =  2t(l)  -  X,  Q  =  h^{l)  -  X,  i|f(l)  =  -E  (the  Euler-Mascheroni 
constant  ,5772)  and  again  X  =  In  sin  0/2.  It  should  be  noted  that 
the  series  S  , S  , S   , S  ,  and  S  _  contained  in  the  scattering  amplitude 

* 
are  not  present  in  the  cross  section,  for  when  u^u,   is  calculated 

these  series  appear  in  terms  of  the  form 


(72)   Tia^  s(M(l)e^''^  -M(1)  e"^*"^]   2A^S 


(X-2^l^(l] 


Since  this  calculation  is  exact  in  the  limit  a  ->  0,  -  —  1  consistency- 
requires  that  such  terms  be  dropped  unless  the  wave  function  (cross 

5  6 
section)  is  known  to  order  cc  (a  )  . 

We  have  calculated  the  cross  section  using  these  formulas 

(equations  (71a-d))  for  scattering  by  lead  (Z  =  82,  a  =  .598)  at  a 

bombarding  energy  of  10  Mev  (p=.9976).   The  results  are  plotted  in 

Figure  3  along  with  the  exact  results  of  Doggett  and  Spencer. 

k.      The  general  form  of  the  coeff icent  of  a  . 

In  general  we  can  write  that 


00   00    r  S 
(75)  \J{9,a,y)   =      2ik  u  (e,a,p)  =  XI   H   frgT   U^  3(6,0,0) 

r=o  S=o    ■  *     •* 


when 

r+S 

(7M  u  .  (e,o,o)  =^-r^  (e,a,v) 

'  hex  b-y 

^  "  ■y=0,a^O 


With  an  increasing  amount  of  analytic  labor  one  can  find  the  U 

r,  b 

by  expansion  of  the  phase  shifts  as  sums  of  the  form: 
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RATIO  OF  RELATIVISTIC     CROSS -SECTION 

TO  RUTHERFORD    SCATTERING   CROSS  -  SECTION 


z  =  82,  energy  =  lOmev 
•    Doggett  a  Spencer 


'A  =  difference 
between  two  curves  above 


0              20*         40*        60*         80*         100*        120*         140*        160*        180* 
9  .- 


Figure  3 
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00 

y     f(n)(P  (cos  ©)  +  P  ^(cosG)) 

^-—^  n       —  n-1 

n=l 


This,  of  course,  was  the  starting  point  of  the  work  descrihed  in  the 
previous  sections.   By  means  of  contour  integration  these  s\jms  were 
expressed  in  the  form  of  integrals  some  of  which  can  be  evaluated  in 
closed  form.   We  shall  now  consider  an  alternative  approach  which  by- 
passes the  intermediate  step,  namely,  expansion  of  the  phase  shift  in 
powers  of  a  and  -^ .      The  total  scattering  amplitude  will  first  be  ex- 
pressed as  a  contoiir  integral  in  the  complex  v-plane;  integral  represent- 
ations will  be  used  for  certain  functions  and  then  derivatives  with 

respect  to  a  and  y   will  be  taken.   The  result  will  be  that  the  U   (©,0,0) 

r,  b 

can  be  expressed  in  terms  of  two  dimensional  real  integrals.   By  a 

simple  change  in  variables  these  integrals  can  be  as  integrals  over  a 

square  in  a  two-dimensional  space. 

k.l     Uon-relativistic  scattering  by  an  inverse  cube  lav  force. 

Before  considering  the  relativistic  Coulomb  problem  itself  we 
shall  digress  briefly  to  discuss  a  related  problem   It  is  well  known 
that  the  differential  equation  describing  the  motion  of  a  relativistic 
classical  particle  moving  in  a  Covilomb  field  differ  from  that  of  a 
non-relativistic  particle  in  that  the  first  contains  an  apparent 
force  term  proportional  to  l/r  .   It  would  appear  reasonable  to  expect 
then  that  there  will  be  mathematical  similarities  in  the  scattering 
amplitude  for  a  particle  obeying  the  non-relativistic  Schrodinger 
equation  moving  in  an  inverse-cube-law  force  field  and  the  scattering 
amplitude  dealt  with  in  this  paper. 

Formally  one  has  the  following  differential  equation: 
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(75) 


r 


2m 


with  the  boundary  conditions: 


X  finite  at  the  origin  , 


(76) 


X   ^>  e^^"  +  h(e) 

r->oo 


ikr 


The  solution  is  given  by 


[^5] 


(77) 


•i5 


X     =     Y.   (2n+l)  P^(cose)  e   ""  j   (kr) 


where 


\--2 


kr^f-     •?  -    (n^y] 


pj—  J      ,  (spherical  Bessel  function) 
n    ' 


T^   =8jt  my'/  h   <  -r  (condition  of  regularity  at  origin) 


o^   =   V(n+^)   -  ^   -  o 


and 


1     ^  r"^^\    n 

(78)      ^(®'")=2lk  ^  (2n+l)  P^(cose)  [e      -  ^| 


OD  00   /  .  sr  [— 

(79)    H(e,  n  )  =  2ikh(e,  T,  )  =  XI  (2n+l)P^(cose)  ^  ^ 

n=o  r=l 


, /,   1^2    2  ,   Iv 
\/(n+^)  -n  -(n+^) 


^     (2v+l)P  (-cose)  ^,,12   2 


2  y      sin  TTV 
D 


r=l 


r/   r;2    2    /^   (+iTr) 


dv 


1+0  - 


The  contoTor  D  is  shown  in  Figure  k   where  the  dashed  portions  indicate 
that  the  path  lies  on  the  second  Riemann  sheet  of  the  functions 
\/(v+l/2)  -  T]       .      The  branch  line'  extends  from  v  =  -l/2  -  n  to 
V  =  -  1/2  +  n  and  is  indicated  by  a  double  line. 


Figvi't.  U 


In  the  top  sheet  Im  V(v+l/2)^-  n^  >  0  at  A  and  Im  V(i'+l/2)^-  1^     <  0 
at  B.   The  portion  of  the  path  lying  on  the  top  Riemann  sheet  is  now 
distorted  so  that  it  runs  parallel  to  the  imaginary  v  axis  except  for 
a  small  semi-circle  centered  about  and  lying  to  the  right  of  the  point 
V  =  -  1/2.  We  have  then 


(80)    H(e,^)  =^^ 


+00  +00 


itT 


-00        -00 


\/cosht-cos0 


f     (2v+l)P  (-cose)      00        ,.    o 


where  K(t,ti)    = 


V^2.    2 


'T     +T1 


r=l 
for  T  >  0 


n/t      -Hr,^     +  (-T^)      for  T  <  0. 


[k(t,ti]' 


dtdT 


[/(v+l/2)2   -^^-{v+l/2)^ 


dv, 


-  kl 


The  contour  A  is  the  integral  which  lies  in  the  bottom  sheet .  Upon 
inspection  we  see  that  a  formal  expansion  of  this  integral  in  powers 
of  T|  yields  terms  of  the  form 


(81)    T1 


2n 


(2V-H)   ^v^-^°^Q) 
"^"  '^^   (v+  1/2)^ 


dv  =  0 


2   1 
as  no  poles  of  the  integrand  are  inclosed  by  the  contour  (for  (3  <  r) ■ 

2   1       2   1 
Note  the  discontinuous  changes  that  occur  for  r\     =  j-     and  r\     >  r  j 

values  for  which  the  scattering  is  -undefined. 
Consider  now  the  function 
oo 


5(t,ti)  =  ^   -7 

r=l 


v^ 


T   +  ri 


We  desire,  in  general,  the  m   drivative  of  this  function  with  respect 

to  n   evaluated  at  n  =  0.   Using  the  integral  representation*-  -^  . 
^ ,00 

(82) 


^ 


T  =  n 


■tX 


J^(tiX)  d\ 


we  find 


(85a)     ^  g(^^^) 
N  2m 

Of) 


n=o 


m 

1 

rq?l 


m     r 
^•'  K      ^    ^^r     D(2m-r,r) 


(83b) 


s2m+l 


2m+l 


g(t,i 


=  0 


where 
(84) 


D(2m-r,r)  =  J2 


P^+Po- • .+P^=2m-r 


C(2m-r;p^,P2, ...p^)  , 


p.  odd 
1 


k2   - 


(85)  C(2m-r,p^...p^)  = 

and 

(P,) 

(86)  J^  "  (0)  = 


(PJ     (Pp)      (Pj 
J^  ^  (0)  J^   (0)...J^  ^  (0) 


Pi  Pg. 


-  ^1^^) 


dz 


Despite  their  seeming  complexity  the  D(2m-r,r)  are  simply  rational 
fractions .  For  t  <  0  we  woiild  find  a  similar  result  for  the 
derivatives  of  -   /  t"  +  ti  +(-t)  except  that  t   in  85a 

would  be  replaced  by  -(-t)  which  is  just  equal  to  t  .   Thus 


(87) 


_1_  cJ^"H(e,^)       ^   \/2  i 

;^  2m  2Tt 

On       T\=o 


+00    +OD 


m; 


(-^) 


r=l 


T 


2m-r-l 


X  D(2m-r,r) 


itT 


/cosht-cos0 


dtdT 


For  m  >  1,   we  have,,    using 

(         itT 
(88) 


2rt  2m-r-l 


(2m-r-2).'(it) 
0 


,,.^.2m-r-2 


t  >  0 
t  <  0 


that 


(89)     ^,     ^ 


2m 


Sn!      N   2m 
On 


H(e,n) 


.A5     V-     /-      xr   (2m-r-2).'      ^,_  w  .  x2m-r-2 

V2     2_     (-^)     ^ „.    '^        D(2m-r,r)(i) 

r=l 


X 


2m-r-2 


\/cosht-cose 


dt 


-  k^ 


and  for  m  =  1  that 


(90) 


i_    o    H(e,Ti) 


2: 


bv' 


-jtiD(l,l) 
CSC  e/2 


(1)  2 

Since  D(l,l)  =  J        (O)   =  l/2  the  cross  section  to  first  order  in  n 

would  be  given  by 


(91) 


kk' 


<^(9)  = 


2  2 

7t   Tl 


1+C 


"e/2 


h-  .2   The  general  terms  in  Mott's  expansion. 

In  treating  the  relativistic  Coulomb  scattering  amplitude  a  similar 
procedure  can  be  carried  out;  the  exact  form  of  the  answer  can  be  in 
several  equivalent  forms  depending  on  the  manner  one  uses  to  separate 
the  phase  shifts  into  distinct  terms .   We  shall  discuss  the  f orm- 
iilation  used  by  Mott  (cf .  equations  (l^)  and  (15))  for  two  reasons: 
1.)   it  can  be  more  easily  compared  with  previous  work,  and  2.)  the 
notation  is  more  compact. 

The  functions  F(9)  and  G(0)  were  defined  in  equations  (l4a,b,c). 
These  may  be  expressed  as 


(92a) 


F(e)  = -i  y     (-1)^  nC  |P  (cose),  -  P  ,(cose)1 
^  '   2   ^j         ^  L_^      l    n-1     _J 


(92b) 


00  2 

G(e)  =  I  Z  (-i)""  n  C^p^(cose)  +  p^_^(coseyj 

n=l  '—  —^ 


Denoting  by  C   the  function  -e""""^"^  -^ — '^ — ~   '^'-     and  setting 

Rv+i+ir) 


-  kh 


(93)  D    =  c     -  c° 

n         n         n 


one   can  rewrite   (92)   as 

(94a)       F(e)   =F     +F,    =F     +-^1"  (-1)^     nD    |p  (cosG)     -   P     ,(cose)1 
'  o         1         o       2  ■^—  ^      '  n  [_n^  '  n-V         _^ 

00  p 

(94b)        G(e)   =G     +G,    =G     +i     y     (-1)°  n     D    [P  (cosG)    -   P     .(cose)"] 
o  1  o       2     ^—^    ^      '  n     n^  '  n-1^  ' 

n=l  '—  — ' 


where 


(95) 


and 


O  2k  |— 1/^       .       ^ 

I  (1-17) 


(96)  G     =  i7F     ctn^9/2    . 

^  00' 

1     ^ 
It   should  be  noted  here  that  Mott  drops  the   siitn  -  ^-r—  ^  (2n+l)P  (cos0) 

n=o 

which  is  zero  if  0  differs  from  zero.   Mott  obtains  the  values  of 

F  and  G  by  physical  considerations,  namely  by  demanding  that,  in  the 

limit  c  ->  00 ,  a  ->  0,  7  finite,  the  scattering  amplitude  reduces  to 

the  non-relativistic  value.   Inasmuch  as  C  — r — ■>  —  and  as  thus 

n  n->QD   n 

2 

n  C   — r >  -n,  this  sum  must  be  retained  in  a  formal  evaluation 

n  n->aD 

of  G  . 

o 

This  point  requires  further  discussion,  for  similar  considerations 
arise  in  more  complicated  fashion  whenever  expansion  of  the  wave 
function  are  made  by  means  of  the  Watson  transformation  in  terms  of 


1.5  - 


00 

Legendre  functions  of  non-integral  order.   The  sum  ^  (2n+l)  P  (cos0) 

n=o 

represents  symbolically  the  angular  dependence  of  the  outgoing  wave 

portion  of  the  incident  wave.   Thus  one  can  find  in  some  texts  the  ex- 

■   [^5] 
pression '-    -' 


(97)  e^^-  ^°^^  =  ^  (2n.l)i°  jjkr)   P^   (cose)     ^::^>  H     ffg  pJcosG) 


X 


ikr        /    -,  vn  -ikr 
e  -   (-1)   e 


This   can  he  written  as 

(98)      e^^^  "°"®     -^— >     ^    le^^''  y  (2n+l)P   (cose)P  (cosO) 
^^    '  r->oo        2ikr    l_         •^—  ^  '   n^  '   n' 


-ikr 


-e"         y  (2n+l)P   (cose)P  (cosir) 


1  ikr   _  /         ^  _ \ 

=  -r^ —     e  5(cos9  -   cos   0) 

ikr 


■e  5(cose  -    cos   rt) 


This  is  clearly  incorrect;  the  basic  error  comes  from  using 

(99)    J^(kr)  ^-^^>     ^       cos  (kr  -  (n+l)rt/2). 

However,  no  matter  how  large  kr  is  made,  there  are  values  of  n 
approximately  equal  to  kr.   It  is  for  these  values,  which  sum  to 
^ikr  cose^  ^^^  ^^^   asymptotic  expression  (99)  is  incorrect.   In 
treating  the  scattering  amplitude,  however,  such  an  expression  can 
usually  be  used  because  the  principal  contribution  to  the  difference 
between  the  total  and  incident  waves  comes  from  partial  waves  of  value 


he 


n(n+l)  ^   kb 
where  b  is  the  wave  number  and  b  the  radius  of  the  region  in  which 
the  potential  energy  is  significantly  comparable  to  the  total  energy. 
In  other  words,  in  non-relativistic  quantum  scattering. 


Yi   (2n+l)  P  (cosG)e  "^ 


2ikr 


will  not  sum  to  the  asymptotic  form  of  the  total  wave  but 


ikr 
e 

2  ikr 


YL   (2n+l)P  (cose)P^^  -  l] 


will  give  the  outgoing  scattered  wave . 

We  shall  now  proceed  to  find  the  formal  expression  of  the  functions 
F  and  G  in  powers  of  a  and  p  =  v/c  .   For  convenience  we  consider 
a  and  y  =   a/p  to  be  independent  variables  and  the  formal  limit  a  -> 

[ke\ 

not  to  imply  7  ->  0.  First  we  express  P  as '-  -" 

(100)  p  =  Lse  +  ^i£^  §3^    p  , 

n     I       n   del      n-1 

This  is  done  in  order  that  we  can  use  the  same  integral  representation 
employed  previously  in  this  paper  despite  a  change  in  contoiir.   Then 

(101)  F  =  +  ^      I  ^      3,^  ^(,,1)    [Tl-cose)v  -  Sine  3gp^.^(-co.e)D(v)dv 

I 

where  C  Is  the  integral  shown  in  Figure  5  and  again 


hi 


Figure  5 


again  the  dotted  contour  P  indicates  a  path  in  the  second 

/~2    2 
Riemann  sheet  of  the  fiinction  V^   ~  ct- 

t 
Distorting  the  contour  C  so  that  it  lies  parallel  to  the 

imaginary  v  axis,  we  find^  employing  (83),  that 

+00 


(102) 


1    _!_ 

2m.'   N  2m 
oa 


2  i  T—  1_  D(2m-r,r) 


X 


|7l-cose)(l/2  +  ix)  -  sine^J      ^itT 


(t-  i/2) 


2m-r 


\/cosh  t  -  COS! 


X  {^-lif{l-e-^f^''     e^^^'^-^^-^^^^d^dtir 


and  that 


-  k8 


(103) 


1       1       S 


2m,  c      F, 


2ml       il      ^   2m       'Tl 
da  ay 


sj2_i_    y     gV2m-r,r)        ,    . 


cx=o,7=o 


ri    s.'(£-s) 


+(X)        +00     OD 


X 


[l-cose)(l/2  +Tljt)   -   sine  ^      .^^^_ 


0-5/2 


-OD        -CX3     o 


(T-    i/2) 


2m-r 


X      [ln2sinhFl    "    [rt-i{]   d^dtdx    . 


Here 


(101.) 


E(s) 


dx^     Rl+x) 


and  E(o)  =  1,  E(l)  =  -^(l)   =  E  (the  Euler  constant  .577...).  The  use 
of  the  following  recursion  relation  is  useful  to  compute  E(s)  for 
large  values  of  s  : 


(105)  E(s+l) 


dx^     r(i+^) 


4'(l+x) 


.  -  E    (p)  '=(-'='Vi'^' 

p=o 


The  ii     -,(1)  s  are  the  polygamma  functions, 
p+i 


In  obtaining  these  last  two  results  we  have  used 


M 


00 

(106)   r(v-ly)    ^_1 B(v-i7,l+2i7)  =  i I  e-^(^-^^V-e"^)^'''d^ 

r(v+i+ir)     rfi+2i7)  Pf 1+217  1  I 


r(i+2i7) 


Re  V  >  0 
where  B  denotes  the  beta  function. 


-  k9 


We  are  now  in  a  position  to  evaluate  the  integral  over  t,    i.e., 
to  sum  the  partial  waves .   The  integrals  to  be  performed  are  of  the 
form 

iT(t-0 


+00 


(107) 


dT 


(T-i/2)P 
and  the  result  of  this  integration  we  shall  denote  symbolically  as 

(108)  2^e-^/2^^-^)iP  5p(t-0 

where,  if  5(t-0  is  the  Dirac  delta  function  and  p  is  positive: 

t-l        u^   U-^ 

(109)  5p(t-0  =   J    j    )    5(u^)du^du2...dUp 

-00  -00 

p  -  fold 


^(t-uP-     H(t-uJ^^'"^'"' 


0 

If  p   is  negative,    the   corresponding  derivative  of  the  delta -function 
is   indicated. 
Finally  then. 


(110) 


h^'h. 


2m- 


\[21  2 


•  2     v^^     V—       D(2m-r)^J:(s)i 


2m- r 


-  :_  rlsU    £-sjI 

r=l     6=0 


+   oo      oo     I — 


X 


(l-cos9)62^_^_^(t-|)-6^_^(t-5)sln©^ 


(.n-ii: 


X 


£n2sinh  ^ 


-s  -t/2 

e 


d^dt 


^  \ycosht-cos( 


and 
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(ill) 


1  ^  ^1 


2ml     n       ^2m^   i 
da     07 


\/21  2 


a?=0,  7=0 


D(2m-r)E(s)l 
tIsI{  P,  -s 


2m- r 


X 


(^^-o.9)5^_^JtH)   +  6^_^_^(t-e)sin0  ^ 


(«-i|)- 


[ 


>C      I     to2sinh  -^ 


■s  e     ' 


\/cosht+cosO 


d|dt   . 


In  general  the  form  of  the  integrals  vlll  be; 


00   t      /      k       k  k 

e^^^(t-l)  ^(iTt-l)  ^(  M2slnh|/2)  ^ 


(cosht-cosO 


J72" 


d|dt  ; 


(those  involving  the  delta  function  and  its  derivative  vill  be  simpler^ 
The  substitutions  y  =  e   ,  x  =  e   brings  this  into  the  form 


n   .^"^r   n'=2r   ,  ,7,^ 


1   1 


(l:t) 


?n  ^ 


«n  y 


1-y 


\/(x-cose)^+  sin^O 


^  dx 


As  previously  noted,  the  dilogarithm  ff^     is  defined  as 


(112) 


^o(z) 


^n(l-0  ^ 
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A  generalization,  called  the  n-logarithm,  is  defined  by  ^  J 


n(^'  -    *    I   ^„-l(S'  -I      ■ 


We  may  further  generalize  these  as  follows: 


^2,i(^)  =-  J    Dnd-lfl'^^  ) 


^n,.(^)=      ^„-l,.(^)1 


Integration  by  parts  allows  us  to  express  oC  ,  (z)  as 

n,k 


^n..^^)=   ^-^n-l,.(^)  -  #^  ^n-2,.---(-)^"C^^2,J^) 


Finally  then,  the  angular  dependence  of  the  correction  terms  to  the  scattered 
wave  will  be  given  by  integrals  of  the  form 

.^  [foxj^CMi-xO  %  vd-^) 

(113)  I     -^^ dx     . 

/  2  2 

V  (x-cosO)    +sin  0 
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5-   Conclusions. 

In  this  paper  we  have  used  the  techniques  of  contour  integration 
to  obtain  in  closed  form  the  differential  cross  section  for  relativistlc 
Coulomb  scattering  up  to  the  fifth  order  in  the  fine-structure-constant. 
The  functional  form  of  the  relativistic  corrections  to  the  scattering 
amplitude  corresponding  to  an  arbitrary  order  of  the  fine  structure 
constant  were  found  in  terms  of  two  dimensional  integrals  involving 
elementary  transendental  functions. 
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Appendix  I 

Evaluation  of  the  Sums 

In  the  evaluation  of  the  sums  listed  in  'teble  1  integrals  of 
the  form: 

^    ^^^^  dt  ,   f(t)  real,  (A  ) 

\/cosh  t  -  cos  0 
o   ^ 

arise.   The  substitution  x  =  e   casts  the  integrals  into  the  form: 

1 

f (-in  x)  dx  (j,\ 

/        2    2 
)       \/(x-cos0)  +sin  0 
o   ^ 

Inasmuch  as  the  integrals  (A  )  are  real  one  must  take  care  that  no  imaginary 

terms  appear  in  the  final  result.   In  particular  the  logarithm  occurring 

t 
in  the  integrals  (B  )  must  be  considered  to  be  in|x|.   The  dilogarithm 

and  n-logarithm  functions  mentioned  in  the  text  require  similar  treatment. 

The  following  functional  relations  were  used  in  simplifying  the  final 

results : 

2 
/^(x)  +  ^f^Cl-x)  =  f-  -   £n|x|in|l-x|  . 

These  relations  differ  from  those  usually  quoted  for  the  reasons  stated. 

The  derivatives  of  the  Legendre  functions  with  respect  to  order  were 
calculated  by  differentiating  the  appropriate  hypergeometric  series  term 
by  term. 

For  example. 


-  5^ 


^  Pv(-cose) 


5 /        l+cosOv 

— 2  F(-v,v+l,l, — g ) 


^v' 


1.  Bimrv  ^   r(^+v+l)r(r-v)  (,^,2^/2)^ 


r=l      (r.-)^ 


=  -2 


r=l    (r!)      ^  —J 


.  _2  f;   (eos^e/2)^     .  J-  ,   2 


r=l      r 


2      =  -2  5C2(cos^e/2) 


Here  we  have  used  the  functional  relation: 

t^(r+l)  =  t-^(r)  +  i  . 

The  corresponding  relationship  for  the  polygamma  functions  follow 
by  differentiation  with  respect  to  r. 
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